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Abstract. Using the new data from the OECD-WTO world network of economic activities we construct 
the Google matrix G of this directed network and perform its detailed analysis. The network contains 58 
countries and 37 activity sectors for years 1995 and 2008. The construction of G, based on Markov chain 
transitions, treats all countries on equal democratic grounds while the contribution of activity sectors is 
proportional to their exchange monetary volume. The Google matrix analysis allows to obtain reliable 
ranking of countries and activity sectors and to determine the sensitivity of CheiRank-PageRank commer¬ 
cial balance of countries in respect to price variations and labor cost in various countries. We demonstrate 
that the developed approach takes into account multiplicity of network links with economy interactions 
between countries and activity sectors thus being more efficient compared to the usual export-import anal¬ 
ysis. The spectrum and eigenstates of G are also analyzed being related to specific activity communities 
of countries. 

PACS. 89.75.Fb Structures and organization in complex systems - 89.65.Gh Econophysics - 89.75.He 
Networks and genealogical trees - 89.20.Hh World Wide Web, Internet 


1 Introduction 

The recent reports of the Organisation for Economic Co¬ 
operation and Development (OECD) [1] and of the World 
Trade Organization (WTO) [2] demonstrate all the com¬ 
plexity of global mannfactoring activities, exchange and 
trade in the modern world. This complexity is rapidly 
growing with time and now it becomes clear that tra¬ 
ditional statistics are increasingly unable to provide all 
the necessary information. Applying modern mathemati¬ 
cal tools and methods to new data sets can allow to under¬ 
stand the hidden trends of the world economic activities. 
Thus the matrix tools for analysis of Input-Out trans¬ 
actions are broadly used in economy starting from the 
fundamental works of Leontief |3l4j with their more re¬ 
cent developments described in [5] . In the last decade the 
development of modern society generated enormous com¬ 
munication and social networks including the World Wide 
Web (WWW), Wikipedia, Twitter and other directed net¬ 
works (see e.g. [6]). It has been found that the concept of 
Markov chains provides a very useful and powerful math¬ 
ematical approach for analysis of such networks. Thus the 
PageRank algorithm, developed by Brin and Page in 1998 
[ 7 ] for the WWW information retrieval, became at the 
mathematical foundation of the Google search engine (see 
e.g. [8]). This algorithm constructs the Google matrix G 
of Markov chain transitions between network nodes and 
allows to rank billions of web pages of the WWW. The 


spectral and other properties of the Google matrix are 
analyzed in [9]. The historical overviews of the develop¬ 
ment of Google matrix methods and their links with the 
works of Leontief are given in [TqITT] . 

The obtained results demonstrate the efficiency of the 
Google matrix analysis not only for the WWW but also 
for various types of directed networks [9]. One of such 
examples is the World Trade Network (WTN) with multi¬ 
product exchange between the world countries. The data 
of trade flows are available at the United Nations (UN) 
COMTRADE database [12] for more than 50 years. The 
results presented in |13ll4j for the WTN show that the 
Google matrix analysis is well adapted to the ranking of 
world countries and trade products and to determination 
of the sensitivity of trade to price variations of various 
products. The new element of such an approach is a demo¬ 
cratic treatment of world countries independently of their 
richness being different from the usual Import and Export 
ranking. At the same time the contributions of various 
products are considered being proportional to their trade 
volume contribution in the exchange flows. 

Here we use the Google matrix analysis developed for 
the multiproduct WTN [14] showing that it can be di¬ 
rectly used for the World Network of Economic Activi¬ 
ties (WNEA) constructed from the OECD-WTO trade in 
value-added database. In a certain sense activities (or sec¬ 
tors) are correlated to products in the WTN. However, for 
the WTN there is exchange between countries but there 
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is no exchange between industries and commodities. Thus 
in [14] it was argued that certain economical features are 
not captured by the COMTRADE database since in real 
economy the traders are industries, not countries; in par¬ 
ticular certain products are transferred to each other (e.g. 
metal and plastic are used for production of cars). In con¬ 
trast to that, the OECD-WTO WNEA incorporates the 
transitions between activity sectors thus representing the 
economic reality of world activities in a more correct man¬ 
ner. 



Fig. 1. World map of countries with color showing country 
import (top panel) and export (bottom panel) with economic 
activity (trade) volume expressed in billions of USD and given 
by numbers at color bars; the gray color marks countries at¬ 
tributed to the ROW group (rest of the world) with exchange 
values 733 (Import) and 1018 (Export) in billions of USD. The 
data are shown for year 2008 with Ac = 57 -h 1 countries (with 
ROW) for the economic activities in all Ns = 37 sectors. Coun¬ 
try names can be found in Table and in the world map of 
countries [22]. 


We note that there has been a number of other inves¬ 
tigations of the WTN reported in |15I16I17I18I19I2Q1^ . 
However, in this work we have the new important ele¬ 
ments, introduced in mni: the analysis of PageRank 
and CheiRank probabilities corresponding to direct and 
inverted network flows and related to Import and Export; 
democratic treatment of countries combined with the con¬ 
tributions of sectors (or products) being proportional to 
their commercial exchange fractions. We point that the 
OECD-WTO TiVA database of economic activities be¬ 
tween world countries and activity sectors has been cre¬ 
ated very recently (2013) and thus this work represents 
the first Google matrix analysis of these data. We stress 


that the usual Import-Export ranking of commercial flows, 
shown in Fig. is not able to take into account all the 
complexity of chains of links between various countries 
and various activity sectors. In contrast to that the ap¬ 
proach developed here takes all of them into account due 
to the powerful method based on the Google matrix. 


2 Methods and data description 

Here we describe the data available for the OECD-WTO 
TiVA network and the mathematical methods used for the 
analysis of this network. The list of Nc = 58 countries (57 
plus 1 for the Rest Of the World ROW) is given in Table^ 
with their flags. Following m we use for countries ISO 
3166-1 alpha-3 code available at Wikipedia. The list of 
sectors with their names is given in Table [^. The fractions 
of sectors in the exchange volume are given in Table for 
years 1995, 2008. 


2.1 Google matrix construction for the OECD-WTO 
WNEA 

We use the OECD-WTO TiVA database released in May 
2013 which covers years 1995, 2000, 2005, 2008, 2009 with 
the main emphasis for years 1995and2008 (2009 data are 
affected by the global crisis and may not be represen¬ 
tative). The network considers Nc = 58 world countries 
given in Table In fact, there are 57 countries and the 
rest of the world, which includes the remaining countries of 
the world forming one group called ROW. There are also 
Ns = 37 sectors of economic activities given in Table 
The sectors are classified according to the International 
Standard Industrial Classification of All Economic Activ¬ 
ities (ISIC) Rev.3 [23]. Here we present results for all 37 
sectors of Table noting that the sectors s = 1,2, ..20 
represent production activities while s = 21,...,37 repre¬ 
sent service activities. The transect ions between service 
sectors are hard to exctract and the future improvements 
of this part of TiVA database are desirable. 

For a given year, the TiVA data extend OECD In¬ 
put/Out tables of economic activity expressed in terms of 
USD for a given year. From these data we construct the 
matrix Med ss' of money transfer between nodes expressed 
in USD: 

Mcd,ss' = transfer from country c', sector sMo c, s (1) 

Here the country indexes are c, c' = 1,. .., Ac and activ¬ 
ity sector indexes are s, s' = 1,..., A^ with Ac = 58 and 
Ns = 37. The whole matrix size is A = Ac x A^ = 2146. 
Here each node represents a pair of country and activity 
sector, a link gives a transfer from a sector of one country 
to another sector of another country. We construct the 
matrix Mcd,ss' from the TiVA Input/Output tables us¬ 
ing the transposed representation so that the volume of 
products or sectors flows in a column from line to line. In 
the construction of Mcd,ss' we exclude exchanges inside 
a given country in order to highlight the trade exchange 
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flows between countries (elements inside country are ze¬ 
ros). 

The ISIC Rev.3 classification of sectors have a sig¬ 
nificant correlation with the UN Standard International 
Trade Classification (SITC) Rev. 1 of products used in 
[14]. There is a clear relationship on the production side 
between ISIC sectors and products of the world exports 
(but not at import level: if all agricultural exports are pro¬ 
duced by the agricultural sector, agricultural products will 
be imported by manufacturing industries such as food pro¬ 
cessing of textile and clothing). There is also another im¬ 
portant difference: the transfer matrix from COMTRADE 
is diagonal in products [T4| (thus there is no transfer from 
product to product), while for the TiVA data there are 
transitions from one sector to another sector and thus the 
matrix of nominal values, in current prices, 0 is not di¬ 
agonal in s, s'. 

For convenience of future notations we also define the 
value of imports Vcs and exports f/* for a given country 
c and sector s as 

Ves = = E ( 2 ) 

c',s' c',s' 

The import Vc = Vcs and export V*g values 

for countries c are shown on the world map of countries 
in Fig. 12 for year 2008. We note that often one uses the 
notion of volume of export or import (see. e.g. El) but 
from the economic view point it more correct to speak 
about value of export or import. 

In order to compare later with the PageRank and 
CheiRank probabilities we define exchange value ranks in 
the whole matrix space of dimension N = Nc x Ng. Thus 
the ImportRank (P) and ExportRank (P*) probabilities 
are given by the normalized import and export values 

Pi = V,,lV, P* = v:JV, ( 3 ) 

where i = s {c — l)A^s, i = and the total 

exchange value is V" = J2c,c',s,s' ^cc',ss' = Ec,s = 

Y V* 

A^cs cs’ 

The Google matrices G and G* are defined as N x N 
real matrices with non-negative elements: 

Gij = aSij + {l-a)viej , G\j = aP*,, + (l-a)<e,-, (4) 

where N = Nc x Ng, a G (0,1] is the damping fac¬ 
tor (0 < (a < 1), Cj is the row vector of unit elements 
{cj = 1), and Vi is a positive column vector called a per¬ 
sonalization vector with Vi = 1 |8ll4j . We note that 
the usual Google matrix corresponds to a personalization 
vector Vi = Ci/N with = 1. In this work, following 
EE], we fix a = 0.5 noting that a variation of a in a 
range (0.5, 0.9) does not significantly affect the probability 
distributions of PageRank and CheiRank vectors IHEI. 
The choice of the personalization vector is specified below. 
Following El we call this approach the Google Personal¬ 
ized Vector Method (GPVM). 


The matrices S and P* are built from money matrices 

■^cc',ss' 






j Mce',ss'/Ve>s' if V'a' 7^0 

\ 1/N if V'.' = 0 

r if yt,, ^ 0 

\ i/AT if y!,, = 0 


(5) 


where c, c' = 1,..., Nc\ 5, s' = 1,..., Ng\ i = s + (c— 1)^^; 
i' = s' (c' — V)Ng] and therefore i, i' = 1,..., V. Here 
Vc's' = ^cs ^cc',ss'- The sum of elements of each column 
of P and P* is normalized to unity and hence the matrices 
G, G*,P, P* belong to the class of Google matrices and 
Markov chains. Thus P, G look at the import perspective 
and P*, G* at the export side of transactions. 

PageRank and CheiRank (P and P*) are the right 
eigenvectors of G and G* matrices respectively at eigen¬ 
value A = 1. The equation for right eigenvectors have the 
form 


V V • ( 6 ) 

j 3 

For the eigenstate at A = 1 we use the notation Pi = 
with the normalization ^Pi = ^iP'^i = 1- 
For other eigenstates we use the normalization = 

The eigenvalues and eigenstates of G, G* 
are obtained by a direct numerical diagonalization using 
the standard numerical packages. 


2.2 PageRank and CheiRank vectors from GPVM 

The components of P^, P*^ are positive. In the WWW con¬ 
text they have a meaning of probabilities to find a random 
surfer on a given WWW node in the limit of large num¬ 
ber of surfer jumps over network links |8|. In the WNEA 
context nodes can be viewed and markets with a random 
trader transitions between them. We will use in the fol¬ 
lowing notation of netwrok nodes. We define the PageR¬ 
ank K and CheiRank P* indexes ordering probabilities P 
and P* in a decreasing order as P{K) > P{K + 1) and 
P*(P) > P*(P* + 1) with P,P* = 1,..., V. 

We note that the pair of PageRank and CheiRank vec¬ 
tors is very natural for economy and trade networks cor¬ 
responding to Import and Export flows. For the directed 
networks the statistical properties of the pair of such rank¬ 
ing vectors have been introduced and studied in [24125113] . 

We compute the reduced PageRank and CheiRank prob¬ 
abilities of countries tracing probabilities over all sectors 
and getting Pc = ^gPcs = ^g P {s {c - l)Ng) and 
P* = Peg = ^g P* (s + (c — l)Ng) with the corre¬ 
sponding Kc and P* indexes. In a similar way we obtain 
the reduced PageRank and CheiRank probabilities for sec¬ 
tors tracing over all countries and getting 
Ps = EcPi^Pi^-^)^s) = Ec Pcs and 
P* = Ec^* ~ l)^s) = ^c^cs with their corre¬ 

sponding sector indexes Kg and P*. A similar procedure 
has been used for the multiproduct WTN data El- 
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In summary we have Kg, K* = 1,..., and Kc, K* = 

1,..., A^c- A similar definition of ranks from import and 
export exchange value can be done in a straightforward 
way via probabilities Pg , P* ^Pc^P^^ Pcs ^ Pqs correspond¬ 
ing indexes Kg , P*, K , • 

To compute the PageRank and CheiRank probabili¬ 
ties from G and G*, keeping a “democratic”, or equal, 
treatment of countries (independently of their richness) 
and at the same time keeping the proportionality of ac¬ 
tivity sectors to their exchange value, we use the Google 
Personalized Vector Method (GPVM) developed in [14] 
with a personalized vector Vi in At the first iteration 
of Google matrix we take into account the relative prod¬ 
uct value per country using the following personalization 
vectors for G and G*: 




V* 

^ n<. 


Vi = 




Vi = 




(7) 


using the definitions (© and the relation i — s-\- {c—l)Ng. 
This personalized vector depends both on sector and coun¬ 
try indexes. As for the multiproduct WTN in m we de¬ 
fine the second iteration vector being proportional to the 
reduced PageRank and CheiRank vectors in sectors, ob¬ 
tained from the GPVM Google matrix of the first itera¬ 
tion: 

= , 8 ) 

In this way we keep democracy in countries but keep con¬ 
tribution of sectors proportional to their exchange value. 
This second iteration personalized vectors are used in the 
following computations and operations with G and G* giv¬ 
ing us the PageRank and CheiRank vectors. This proce¬ 
dure with two iterations forms our GPVM approach. The 
difference between results obtained from the first and sec¬ 
ond iterations is not very large (see Figs. |^, but the 
personalized vector for the second iteration gives a reduc¬ 
tion of fluctuations. In all Figures after Fig.j^we show the 
GPVM results after the second iteration. 

As for the WTN it is convenient to analyze the distri¬ 
bution of nodes on the PageRank-CheiRank plane A"*). 
In addition to two ranking indexes K, if* we use also 
2DRank index K 2 which describes the combined contri¬ 
bution of two ranks as described in [25|. The ranking list 
K 2 {i) is constructed by increasing K ^ K-\-l and increas¬ 
ing 2DRank index K 2 {i) by one if a new entry is present 
in the list of first A"* < K entries of CheiRank, then the 
one unit step is done in iC* and K 2 is increased by one if 
the new entry is present in the list of first K < iC* entries 
of CheiRank. More formally, 2DRank K 2 {i) gives the or¬ 
dering of the sequence of nodes, that appear inside the 
squares [1,1; K = A"* = k; ...] when one runs progres¬ 

sively from k = 1 to N. Additionally, we analyze the dis¬ 
tribution of nodes for reduced indexes (A'c, A"*), (A^, A*). 

The localization properties of eigenstates of G, G* are 
characterized by the inverse participation ration (IPR) 
defined as ^ = (X^i Si This quantity de¬ 

termines an effective number of nodes contributing to a 
formation of a given eigenstate (see details in |9]). 



Fig. 2. Dependence of probabilities of PageRank P(A), 
CheiRank P*(A*), ImportRank P(A) and ExportRank 
P*{K*) on their indexes in logarithmic scale for WNEA (or 
OECD-WTO TiVA network) in 2008 with a = 0.5, Nc = 58, 
Ns = 37, V = Ac X W = 2146. Here the results for the GPVM 
after the first and second iterations are shown for PageRank 
(CheiRank) in red (blue) with dashed and solid curves respec¬ 
tively. Probabilities for ImportRank and ExportRank from 
exchange value are shown by magenta and cyan thin curves 
respectively. 



Fig. 3. Probability distributions of PageRank and CheiRank 
for sectors Ps(Ks), Ps{K*) (left panel) and countries Pc(Kc), 
Pc{Kc) (right panel) in logarithmic scale for WNEA (or 
OECD-WTO TiVA network) from Fig|^ Here the results for 
the first and second GPVM iterations are shown by red (blue) 
curves for PageRank (CheiRank) with dashed and solid curves 
respectively (with a strong overlap of curves). The probabilities 
from the exchange value ranking are shown by thin magenta 
and cyan lines for ImportRank and ExportRank respectively. 


2.3 Correlators of PageRank and CheiRank vectors 

As in previous works |24l25ll3j we consider the correlator 
of PageRank and CheiRank vectors: 

N 

/t = Ar^P(i)P*(i) - 1 . (9) 

i=l 

The typical values of k are given in |9| for various net¬ 
works. 
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For the global PageRank and CheiRank probabilities 
the sector-sector correlator matrix is defined as: 


= iVcF 


P{s + (c - l)N,)P*is' + (c - 1)N,) 


Ee' Pis + (c' - l)Ns) Ec" P*is' + (c" - l)Ns)_ 

( 10 ) 


Then the correlator for a given sector is obtained from 
(10) as; 

( 11 ) 


- f^ss' ^s,s' 5 


where Sg^s' is the Kronecker delta. 

We also use the correlators obtained from the probabil¬ 
ities traced over sectors (Pc = Pgc) and over countries 
{Pg = ^sc) which are defined as 


Ne 


Ns 


t(c) = N,J2 PoK - 1. <s) =N,J2 PsP: - 1 ■ (12) 


C=1 


S = 1 


In the above equations the correlators are com¬ 

puted for PageRank and CheiRank probabilities. We can 
also compute the same correlators using probabilities from 
the exchange value in ImportRank P and ExportRank P* 
defined by (§. 

The obtained results are presented in the next Section 
and at the web site [26l • 


3 Results 
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We apply the GPVM approach to the data sets of OECD- Fig. 4. Country positions on PageRank-CheiRank plane 
WTO TiVA of WNEA and present the obtained results {Kc,K*) obtained for the WNEA by the GPVM analysis (top 
below. panels), ImportRank-ExportRank of exchange value (middle 

panels), and PageRank-CheiRank plane of WTN ranking of 
trade in all commodities from [13] (bottom panels) shown for 
Kc,K*c < 60. Left (right) panels show year 1995 (2008). 

3.1 PageRank and CheiRank probabilities 


The dependence of probabilities of PageRank P{K) and 
CheiRank P*{K'^) vectors on their indexes P, P* are shown 
in Fig. [^for a selected year 2008. The results can be ap¬ 
proximately described by an algebraic dependence P oc 
IjK^^ P* (X 1/P*^ with the fit exponent value P = 
0.385 ± 0.014 for PageRank and /3 = 0.486 ± 0.02 for 
CheiRank for P, P* < 10^. In contrast to WWW and 
Wikipedia networks (see e.g. El) there is no significant 
difference of f3 between two ranks that can be attributed 
to an intrinsic property of economy networks to keep econ¬ 
omy balance of commercial exchange. The probability vari¬ 
ation is reduced for the Google ranking compared to the 
value ranking. This results from a “democratic”, or equal 
grounds ranking of countries used in the Google matrix 
analysis. The obtained data also show that the variation 
of probabilities for 1st and 2nd GPVM iterations are not 
very large that demonstrates the convergence of this ap¬ 
proach. 


3.2 Ranking of countries and sectors 

After tracing the probabilities P(P),P*(P*) over sec¬ 
tors we obtain the distribution of world countries on the 
PageRank-CheiRank plane (Pc,P*) presented in Fig. 
for WNEA in years 1995, 2008. In the same figure we 
present the rank distributions obtained from ImportRank- 
ExportRank probabilities of exchange value and the re¬ 
sults obtained in m for the WTN with all commodities. 
For the GPVM data we see the global features already 
discussed in m- the countries are distributed in a vicin¬ 
ity of diagonal Pc = P* since for each country the size of 
imports is correlated with the size of exports, even if trade 
is never exactly balanced and some countries can sustain 
significant trade surplus or deficit. The top 20 list of top 
P 2 countries recover 13 of 19 countries of G20 major world 
economies (EU is the number 20) thus obtaining 68% of 
the whole list. This is close to the percent obtained in 
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m for trade in all commodities. The Google ranking for 
WNEA and WTN (top and bottom panels in Fig.[^ gives 
different positions for specific countries (e.g. Russia im¬ 
proves its position for WNEA with the opposite trend for 
China) but the global features of distributions of WNEA 
and WTN remain similar corresponding to the same eco¬ 
nomical forces. 


PageRank-CheiRank (GPVM) 


ImportRank-ExportRank 
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• 3C15T16FOD ■ 

• 4C17T19TEX 
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16 C32 CMQ 
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-I- 27C65T67FIN ♦ 37 C95 PVH 

28 C70 REA 
^ 29C71RMQ 
^ 30C72ITS 


Fig. 5. Two dimensional ranking of sectors on the (Kg,Kg) 
plane using the GPVM approach for PageRank and GheiRank 
(left panels) and ImportRank-ExportRank (right panels). Each 
sector is represented by its specific combination of color and 
symbol. The list of all 37 sectors are given in Table Top 
panels show the case for the year 1995 and bottom panels for 
the year 2008. 


After tracing over countries we obtain the PageRank- 
CheiRank plane of activity sectors shown in Fig. We 
see that some sectors are export oriented (e.g. s = 2 
CIOTI 4 Mining at = 1 in 2008) others are import 
oriented (e.g. s = 23 C50T52 World Retail and Trade of 
motors etc. at = 1 in 2008). The ImportRanking gives 
a rather different import leader s = 1 C23 Manufacture of 
coke, refined petroleum products etc. with Kg = I in 2008. 
Thus the Google ranking highlights highly connected net¬ 
work nodes while Import-Export gives preference to high 
value neglecting existing network relations between vari¬ 


ous countries and activity sectors. We can also order sec¬ 
tors by 2DRank index K 2 getting for PageRank-CheiRank 
top sectors s = 25,23,8 at K 2 = 1,2,3 while Import- 
Export gives 5 = 8,11,14 for top K 2 values in 2008 (more 
data are given at [26]). We note that s = 2b corresponds 
to Transport which has many network connections thus 
taking the top K 2 position. We note that asymmetry of 
ranking of products has been discussed in [14] for COM- 
TRADE data, however, the comparison with these data is 
not so simple since the correspondence between products 
and activity sectors is not straightforward. Of course, for 
the WNEA the asymmetry of sector ranking exists even 
for Export-Import ranking, in a drastic difference from 
the WTN, since there are interactions between activity 
sectors. 

The global ranks of top 20 countries and their activities 
are given in Table for 2008. The top 3 places of PageR¬ 
ank K = 1,2,3 are taken by Germany {Manufacture of 
motors etc. s = 18), USA {Public administration and de¬ 
fence s = 33), ROW (also s = 33). Thus imports of arms 
and weapons play a very important role. In contrast for 
ImportRank K = 1,2,3 we find rather different results 
with USA (petroleum s = 7), Japan (also s = 7), and 
only then USA {s = 33). For CheiRank = 1,2,3 we 
find ROW, Russia, Saudi Arabia (5 = 2 CIOTI 4 Mining) 
while for ExportRank we have ROW, Saudi Arabia, Rus¬ 
sia (s = 2 CIOTI 4 Mining) respectively. Thus Russia goes 
ahead of Saudi Arabia due to a broad network of activ¬ 
ity and trade connections (a similar effect has been found 
in |13ll4j for trade in petroleum). The top 3 positions of 
2DRank K 2 = 1, 2, 3 are taken by Germany {s = 8 Man¬ 
ufacture of chemicals etc.). USA {s = 27 Finance etc.), 
Germany {s = 13 Manufacture of machinery etc.). 

We can fix a certain activity sector s and then consider 
local ranking of countries in {Kc, Kc"^) plane. Three exam¬ 
ples are shown in Fig.j^for s = 21 {Electricity, gas, water), 
28 {Real estate activity, 1 {Agriculture). The comparison 
of Google ranking (left column) with value Import-Export 
ranking (right column) shows importance of network con¬ 
nections highlighted by the GPVM, thus Russia moves 
from Kc^ = 4 on right panel to = 2 on left panel 
for s = 21 due to its broad links with Europe and Asia. 
For s = 1 case in bottom panels of Fig. we find that 
the Import-Export ranking distribution is more else to di¬ 
agonal comparing to the PageRank-CheiRank case that 
we attribute to effect of indirect links present in the later 
case. 

The distribution of nodes on the global {K, i^*) plane 
is shown in Fig. for Google ranking (left panel) and 
Import-Export ranking (right panel) in 2008. The major¬ 
ity of countries are shown by gray squares while 6 selected 
countries are marked by colors. The comparison of two 
panels show that in the Google ranking the positions of 
USA are improved (more black symbols at top K 2 posi¬ 
tions) while for China the positions (green symbols) are 
weakened. We attribute this to a broader network con¬ 
nections of USA in important activity sectors world wide 
(e.g. military activities and defense). 
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Fig. 6. Left column panels show results of the GPVM data for 
country positions on PageRank-GheiRank plane of local rank 
values Kc^ Kc^ ordered by (iLcs, K*cs) for specific sectors with 
s = 21 (top), s = 28 (center) and s = 1 (bottom). Right column 
panels show the Import Rank-Export Rank planes respectively 
for comparison. Data are given for year 2008. Each country is 
shown by its own flag as in Eig[^ 


3.3 Correlation properties of PageRank and CheiRank 

The directed networks can be characterized by the corre¬ 
lator n of PageRank and CheiRank vectors. For various 
networks the properties of n are reported in |24l9j . There 
are directed networks with small or even slightly negative 
values of e.g. Linux Kernel or Physical Review citation 
networks, or with ^ 4 for Wikipedia networks and even 
larger values k, ^ 116 for the Twitter network. 

The correlators of WNEA for various sectors are shown 
in Fig. Almost all correlators hZs are positive being dis¬ 
tributed in a range (0,1). A small negative value appears 
only for s = 37 {Private households ete.) corresponding to 
anti-correlation between buyers and sellers. The largest 
correlator Kg is for s = 29 {Renting of maehinery ete.) 
shows that sales of machinery correlates with their pur¬ 
chases probably because components are needed to pro- 



K K 


Fig. 7. Global plane of rank indexes {K,K*) for PageRank- 
GheiRank (left panel) and Import Rank-Export Rank (right 
panel) for N = 2146 nodes in year 2008. Each country and 
sector pair is represented by a gray square. Some countries 
are highlighted in colors : USA in black, South Korea in red, 
Ghina (and Taiwan) in green, Russia in blue, France in yellow 
and Brazil in orange. 



Fig. 8. PageRank-GheiRank correlators Ks from the GPVM 
(see 


( put , (II)) are shown as a function of the sector in¬ 
dex s with the corresponding symbol from Fig[^ PageRank- 
GheiRank and Import Rank-Export Rank correlators are shown 
by solid and dashed lines respectively, where the global corre¬ 
lator n it is shown in black, the correlator for countries k,{c) 
(12) is shown by red lines, the correlators for sectors k,{s) ( p^ 
is shown by blue lines. Here sector index s is counted in order 
of appearance in TableThe data are given for year 2008 with 
Ns = 37, Nc = 58, V = 2146. 


duce machines produced by firms in the same industrial 
sectors. 

The matrix of correlators between sectors s^ s' is shown 
in Fig. for years 1995, 2008. It is interesting to see a 
significant shift of line of maximal correlators located in 
1995 at s' = 28 {Real estate aetivities) to 5 = 29 {Renting 
of maehinery ete.) in 2008. We also see that there are less 
correlations between sectors in 2008 compared to 1995. A 
further more detailed analysis of correlations would bring 
a better understanding of hidden inter-relations between 
various sectors of economic activity. 
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Fig. 9. Top panels: Spectrum of Google matrices G (left) and 
G* (right) represented in the complex plane of A. The data are 
for year 2008 with a = 1, N = 2146, Nc = 58, Ns = 37. Four 
eigenvalues marked by colored circles are used for illustration 
of eigenstates in Fig. |10| and Table Bottom panels: Inverse 
participation ratio (IPR) ^ of all eigenstates of G (left) and 
G* (right) as a function of the real part of the corresponding 
eigenvalue A from the spectrum above. 


3.4 Spectrum and eigenstates of WNEA Google matrix 

The results obtained for the Wikipedia network [29] and 
the multiproduct WTN m demonstrated that the eigen¬ 
vectors of G and G* with large eigenvalue modulus |A| 
select certain specific communities. Thus it is interesting 
to analyze the properties of eigenvalues for the WNEA. At 
a = 1 the gap between A = 1 and other eigenvalues char¬ 
acterize the rate of system relaxation to the equilibrium 
stationary PageRank state (for G). The presence of small 
gap indicates that the mixing and relaxation in the system 
are developed only after many iterations of G matrix (see 
more discussion in 0)- 

The matrix size of WNEA is relatively small and the 
whole spectrum A of G, G* can be determined by direct 
matrix diagonalization. The spectrum is shown in top pan¬ 
els of Fig. It is characterized by a significant gap be¬ 
tween A = 1 and other eigenvalues with |A|<0.7at(a = l. 
We attribute this to a large number of inter-connected 
links between matrix nodes (countries and sectors) which 
is usually responsible for appearance of the spectral gap 
(see m, where the gap increases with the increase of num¬ 
ber of random links per node). We also note that the max¬ 
imal value of \ImX\ < 0.2 is relatively small due to pres¬ 
ence of links going in direct and inverse directions between 
nodes. These features show that the relaxation processes 
to the steady-state PageRank vector are relatively rapid 
on the WNEA. Indeed, the relaxation is governed by the 
exponent exp(—Z\At) where AX ^ 0.25 the gap for for 
WNEA in Fig. [^and t is number of iterations of G. 

The properties of eigenstates are characterized by the 
IPR ^ shown in bottom panels of Fig. [^ We find that the 
main part of states have ^ ^ so that they occupy only 
a small fraction of nodes corresponding to localized states 
(see discussion about the Anderson localization of Google 
matrix eigenstates in |9|28j ). 

The dependence of amplitudes {ipil of a few eigen¬ 
states, ordered by a local rank index Ki corresponding 
to a monotonic amplitude decrease, are shown in Fig. 
The names of top 10 nodes of these eigenstates are given in 
Tabled The red curve in Fig. [Tq| selects mainly the sector 



Fig. 10. Eigenstates amplitudes \'ipi\ ordered by its own de¬ 
creasing amplitude order with local rank index Ki for 4 differ¬ 
ent eigenvalues of Fig. [^(states are normalized as \^i\ = 1). 
The four examples are A = 0.4993 (red), A = 0.3746 + 0.0126z 
(green), A = 0.6256 (blue) and A = —0.0001 + 0.1687i (ma¬ 
genta). Node names (country, sector) for top ten largest am¬ 
plitudes of these eigenvectors are shown in Table 


5 = 4 {Manufacture of textiles etc.) with close links be¬ 
tween China, Italy, USA and ROW; the green one selects 
5 = 18 {Manufacture of motor vehicles etc.) with close 
links between Argentina, Brasil, Japan and Germany; the 
blue state corresponds to 5 = 16 {Manufacture of radio, 
television and communication equipment and apparatus) 
in the Asian region (China, Korea, Chinese Taipei, Sin¬ 
gapore, Malaysia); the magenta state represents sector 
5 = 2 {Mining etc.) with related countries like Russia, 
Saudi Arabia, ROW, Norway. These results coincide with 
the previous observations for Wikipedia-type network [29] 
that the eigenstates of G and G* select specific communi¬ 
ties of the network nodes. Similar properties of eigenstates 
of G of the multiproduct WTN have been found in m- 


3.5 Sensitivity to price variations 

The ranking of WNEA nodes provides interesting and im¬ 
portant information. In addition, the established matrix 
structure of G, G* of WNEA also allows to study the sensi¬ 
tivity of the world economic activities to price variations. 
There are certain parallels with the multiproduct WTN 
analyzed in [14] but there are also new elements specific 
to the WNEA. 

To analyze the sensitivity of price variation in a cer¬ 
tain activity sector 5 we increase from 1 to 1 -b the 
money transfer in the sector 5 in Mccss' where Sg 

is a dimensionless fraction variation of price in this sector. 
After that the matrices G, G* are recomputed in the usual 
way described above and their rank probabilities P, P* are 
determined. Then we compute the derivatives of probabil¬ 
ities of PageRank D = dP/dSg = AP/Sg and CheiRank 
P* = dP""/ddg = AP""/5g. We do these computations at 
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Fig. 11. Left panel: Derivatives D = dPfdSr and D* = 
dP""Iddr for a price variation 67 of 7 C23 PET (Manufae- 
ture of eoke, refined petroleum produets and nuelear fuel) for 
year 2008. Right panel: Logarithmic derivatives Di = D/P 
and D/ = D* /P* for the same case as left panel. Codes in 
panels give sector number s = 1, ...37 described in Table 
country codes are from Table The group of points, high¬ 
lighted by the dashed box, represents 58 nodes of the form 
(country, s = 7) where s = 7 is C23 PET (Manufaeture of 
eoke, refined petroleum produets and nuelear fuel). 


sufficiently small Ss values checking that the variations of 
P, P* are linear in 63 . In addition we also compute the 
logarithmic derivatives Di = dlnP/dSs, D/ = dhiP""/d5s 
which give us relative changes of P, P*. 

The sensitivities to price of s = 7 (Manufacture of 
coke, refined petroleum products and nuclear fuel) are shown 
in Fig. The data for P, P* in the left panel show a 
rather complex picture with a significant derivatives not 
only for s = 7 but also for countries with sectors: s = IS 
(Manufacture of motor vehicles, trailers and semi-trailers) 
at strongly negative P for Germany. USA, Japan; s = 25 
(Land transport; transport via pipelines etc) at significant 
positive P for Germany. Korea, Denmark, Singapore; of 
course, for 5 = 7 we have positive P*, but also for s = 2 
related to mining and negative P* for 5 = 8 (Manufacture 
of chemicals and chemical products) for USA and Ger¬ 
many. The logarithmic derivatives provide strong relative 
changes and are shown in the right panel of Fig. [F 



I_I I 

Fig. 12. Same as the left panel of Fig. |ll| but using proba¬ 
bilities from the trade value. In the right panel, Pf = 0.9348 
if s = s' and Pf = -0.0633 if s / s'. 


A similar analysis can be done using the probabilities 
P, P* from the exchange value probabilities ([^ instead 
of the above PageRank and CheiRank probabilities. The 
results for the value probabilities are presented in Fig. 
for the same case as in Fig. F We see that the results are 
drastically different especially for the logarithmic deriva¬ 
tives Di , D /. In fact Di , D/ cannot give correct picture of 
sensitivity to price variations since for the monetary ex¬ 
change the network links between nodes are not taken into 
account and there is only a mechanical re-computation of 
the value normalization. A similar situation appears also 


and Fig.y^that the Google matrix approach provides new 
elements for the economic activity analysis going signifi¬ 
cantly beyond the usual consideration of Import-Export 
method. 

The new element of the WNEA, compared to the mul- 
tiproduct WTN, is existence of transfers between sectors 
of the same economy. This allows us to consider the sen¬ 
sitivity not only to sectoral prices but also the sensitivity 
to labor cost in a given country c (e.g. price shock af¬ 
fecting all industries in the same country). This can be 
taken into account by the introduction of the dimension¬ 
less labor cost change in a given country c by replacing 
the related monetary flows from coefficient 1 to 1 + ctc in 
Me,,,,,; 0 for a selected country c. 

Of course, the above derivatives over price of activity 
sector and labor country cost give only an approximate 
consideration of effects of price variations which is a very 
complex phenomenon. For an economic discussion of the 
effect of price shocks on international production networks 
we address a reader to the research performed in [30] . We 
will see below that our approach gives results being in 
a good agreement with economic realities thus opening 
complementary possibilities of economic activity analysis 
based on the underlying network relations between coun¬ 
tries and activity sectors which are absent in the usual 
Import-Export consideration. We present the results on 
sensitivity to sector prices and labor cost in next subsec¬ 
tions. 


for the multiproduct WTN m- Thus we see from Fig. 
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3.6 Price shocks and trade balance sensitivity 

On the basis of the obtained WNEA Google matrix we can 
now analyze the trade balance in various activity sectors 
for all world countries. Usually economists consider the 
export and import of a given country as it is shown in 
Fig. El Then the trade balance of a given country c can be 
defined making summation over all sectors: 

Be = 'E{p:,-Pcs)/^{p:s+Pcs) = {p:-Pc)/{p:+Pc). 

s s 

(13) 

In economy. Pc, P* are defined via the probabilities of 
trade value Pcs^Pcs ([^- Iii matrix approach, we 

define Pcs^Pcs PageRank and CheiRank probabilities. 
In contrast to the Import-Export value our approach takes 
into account the multiple network links between nodes. 
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Fig. 13. World map of CheiRank-PageRank balance Be = 
(P* ~ Pc)/{Pc + Pc) determined for all Nc = 58 countries in 
year 2008. Top panel shows the probabilities P and P* given 
by PageRank and CheiRank vectors; the value of ROW group 
is Pc =58 = 0.023. Bottom panel shows the probabilities P and 
P* computed from the Export and Import value; the value of 
ROW group is Pc =58 = 0.16. Names of the countries are given 
in Table and in the world map of countries [22]. 



Fig. 14. Derivative of probabilities balance dBc/dSr over price 
of sector s — 7 C23PET for year 2008. Top panel shows the 
case when Be is determined by CheiRank and PageRank vec¬ 
tors as in the top panel of Fig[^ the value of ROW group 
is dB^s/dSr = 0.04. Bottom panel shows the case when Be 
is computed from the Export-Import value as in the bottom 
panel of Eig[^ the value of ROW group is dB^g/dSr = —0.07. 
Names of the countries can be found in Table [T] and in the 
world map of countries [22]. 


The comparison of CheiRank-PageRank balance with 
Export-Import balance for the world countries shown in 
Fig. for year 2008. Each country is shown by color wh- 
cih is proportional to the country balance Be (13) with 
the color bar given on the figure. For Export-Import bal¬ 
ance we see the dominance of petroleum producing coun¬ 
tries Saudi Arabia, Russia, Norway with the largest val¬ 
ues. The CheiRank-PageRank balance highlights 

new features placing on the top Russia, Norway, Germany, 
China. In fact, USA has now a slightly positive balance 
in top panel of Fig. 13) while it was negative before in 
bottom panel of same figure. We see that the broad net¬ 
work of economic activity relations and links makes the 
economies of the above countries more important in the 
world economy while Saudi Arabia, with the largest pos¬ 
itive Export-Import balance, looses its leading position. 
Indeed, the trade of this country is mainly oriented to 
USA and nearby countries that reduces its importance for 
world economy (a similar effect has been observed with 
COMTRADE data p^TTij b 

The sensitivity of country balance dBc/dd^ to price 
variation of sector s = 7 Manufacture of coke, refined 
petroleum products and nuclear fuel is shown in Fig. 

For Export-Import in bottom panel the most sensitive 
countries are Lithuania (positive) and Vietnam (negative). 
Lithuania does not produce petroleum, but in fact in 2008 


there was a large oil refinery company there which had a 
large exportation value (see e.g. 

http://en.wikipedia.org/wiki/Economy_of_Lithuania). The 
Export-Import approach shows that Russia is slightly pos¬ 
itive, even less positive is Saudi Arabia, China and Ger¬ 
many are close to zero change, USA is only very slightly 
positive. The results of CheiRank-PageRank sensitivity 
(top panel) are significantly different showing strongly pos¬ 
itive sensitivity for Saudi Arabia, Russia and strongly neg¬ 
ative sensitivity for China, Germany and Japan; USA goes 
from slightly positive side in bottom panel to moderate 
negative one in top panel. The CheiRank-PageRank bal¬ 
ance demonstrates much higher sensitivity of Russia, Saudi 
Arabia and China to price variations of s = 7 sector com¬ 
paring to the case of Export-Import value analysis. The 
economies of Germany, China and Japan are also very 
sensitive to petroleum prices that is correctly captured by 
our analysis. We consider that the CheiRank-PageRank 
approach describes the economic reality from a new com¬ 
plementary angle and that provides new useful informa¬ 
tion about complex trade systems. We also note that the 
highly negative sensitivity of China to petroleum prices 
has been also obtained on the basis of Google matrix anal¬ 
ysis of COMTRADE data (see Fig.21 in [T4]L 
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Fig. 15. Derivative of partial probability balance of sector 
s defined as dBcs/dSs' over sector s' = 7 C23PET price Sr 
for year 2008. Here Bcs = (F*. - Pcs)/{Pc + Pc) and s = 2 
(CIOTI 4 MIN, Mining, extraction,...) from Table The sec¬ 
tor balance sensitivity of countries Bcs is determined from 
CheiRank and PageRank vectors (top panel) and from the ex¬ 
change value of Export-Import (bottom panel); the values of 
ROW group are dBr>s, 2 /d 5 r — 0.05 and dBr>s, 2 /d 5 r — —0.03 
respectively. Names of the countries can be found at Table 
and in the world map of countries [22]. 


It is also possible to determine the cross-sensitivity of 
activity sectors to price variation. For that we determine 
the partial exchange balance for a given sector s defined 
as 

Be. = (Pi-Pc.)/5](P*. + a=.) = (P*.-Pc.)/(Pc*+Pc), 

(14) 

SO that the global country balance is Be = Xls Then 
the sensitivity of partial balance of a given sector s in re¬ 
spect to a price variation of a sector s' is given by the 
derivative dBcs/dSg'. The results for s = 2, s' = 7 are 
shown in Fig. We see that two methods give results 
with even opposite signs. According to the Google matrix 
analysis the increase of petroleum prices stimulates devel¬ 
opment of mining while for the Export-Import approach 
the result is the opposite. In our opinion, the absence of 
links and next step relations between countries and sectors 
in the Export-Import methods does not allow to take into 
account all complexity of economy relations. In contrast 
the CheiRank-PageRank approach captures effects of all 
links providing more advanced indications. 

The sensitivities dBc/dSs' of CheiRank-PageRank bal¬ 
ance of China and USA to price variation of sectors s' are 
presented in Fig. We see two rather different profiles. 




Fig. 16. Top (China) and bottom (USA) panels show deriva¬ 
tive dBc/dSg' of country total probability balance Be over 
price Ss' of sector s' for year 2008 (black points connected 
by dashed line); derivatives of balance without diagonal term 
{dBc/dSs' — dBes'/dSs') are represented by open red circles. 
The sector balance of countries Bcs and Be are determined 
from CheiRank and PageRank vectors. The sectors correspond¬ 
ing to sector index s or s' are listed in Tabl^ 


Thus, for China the derivative dBc/dSg' is positive for sec¬ 
tors 5 = 4,14,16 {Manufacture of textiles; office machin¬ 
ery; radio etc,) and negative for s = 7, 25, 27 {Petroleum; 
Land transport etc.; Financial intermediation etc.). For 
USA the sensitivity is significantly positive for s = 23, 29, 32 
{Sale of motor vehicles etc.; Renting of machinery and 
equipment etc.; Other business activities) and negative for 
5 = 11 {Manufacture of basic metals). Thus the economic 
activities of these two countries have very different strong 
and weak points. We note that the sensitivity without the 
diagonal term {dBc/dSs> — dBcs' / dSs ') has negative values 
for almost all sectors for both countries. 

The matrices of cross-sector sensitivity dBcs/dSs' are 
shown for China and USA in Fig. 17 Such matrices pro¬ 
vide a detailed information of interconnections of various 
activity sectors. Thus for USA we see that its 5 = 8 {Man¬ 
ufacture of chemicals etc.) has a significant negative sen¬ 
sitivity to s' = 7, 23, 25 {Petroleum; Renting of machinery 
and equipment etc.; Land transport etc.). Indeed, chemical 
production is linked with petroleum, machinery and trans¬ 
port. For China we find that its sector 5 = 11 {Manufac¬ 
ture of basic metals) has a negative sensitivity to 5 ' = 8, 23 
{Manufacture of chemicals etc.; Renting of machinery and 
equipment etc.)] also s = 14,16 have a negative derivative 
in respect to s' = 11). 

Of course, the cross sensitivity to price variations in 
one sector and their effects on another sector, based on 
(14), is a very delicate thing since a price in one sector 
can affect prices in other sectors also in other manner since 
economic systems learn and adapt while here we consid¬ 
ered only linear algebraic relations without any adapta¬ 
tion features. However, even being linear, the Google ma¬ 
trix approach provides a detailed information on hidden 
interactions and inter-dependencies of various economic 
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Fig. 17. China (top) and USA (bottom) examples of deriva¬ 
tive dBcs/dSg' of partial probability balance Bcs of sector s 
over price Ss' of sector s' for year 2008. Diagonal terms, given 
by y = dBcs/d6s for s = s', are shown on the top panels of each 
example. Sectors s' and s are shown in x-axis and y-axis respec¬ 
tively (indexed as in Tabl^from 1 to 37), while dBcs/dSg' is 
represented by colors with a threshold value given by H-e and 
—e for negative and positive values respectively, also shown 
in red dashed lines on top panels with diagonal terms. Here 
e = 0.01 for USA and China; partial balance Bcs is dehned by 
CheiRank and PageRank probabilities. 



activities for various countries that can provide a useful 
message even for nonlinear adapting systems. 



Fig. 18. Derivative of probabilities balance dBcjdac' over 
labor cost of China c = 37 for year 2008. Top panel shows 
the case when Be is determined by CheiRank and PageRank 
vectors; here the special values are dB^s/dasr = —0.0146 for 
ROW group (gray) and dB^r/da^r = 0.3217 for China (ma¬ 
genta) . Bottom panel shows the case when Be is computed from 
the Export-Import value; the special values are dB^^jda^i — 
—0.0352 for ROW group (gray) and dHsr/dcrsr = 0.4810 for 
China (magenta). Names of the countries can be found in Table 
[^and in the world map of countries [22]. 


3.7 World map of sensitivity to labor cost 

Using the established structure of WNEA we can study 
the sensitivity of country balance dBcjda'^ to the labor 
cost in different countries. At the difference of sectoral 
shocks on one product, here the price shock affects all in¬ 
dustries in a country. As before, the change in price has 
to be small enough for the resulting simulation to remain 
in a neighbourhood of the original data. Indeed, larger 
shocks would trigger a series of substitution effects divert¬ 
ing trade to other partners. 

The derivative dBc/da'^ is computed numerically as de¬ 
scribed in Sec. 3.5. The world sensitivity to the labor cost 
of China is shown in Fig.[^ Of course, the largest deriva¬ 
tive is found for China itself {dBc/dac at c = 37 from 
Table . The effect on other countries is given by non¬ 
diagonal derivatives at c 7 ^ c' = 37. From the CheiRank- 
PageRank balance we find that the most strong negative 
effect (minimal negative dBc/dac') is obtained for USA, 
Germany, UK; a positive derivative is visible only for Chi¬ 
nese Taipei {s = 38) and S.Korea {s = 19). For the 
Export-Import balance the results are rather different: at 
first all derivatives at c 7 ^ c' are negative; among the most 
negative values are such countries as Hong Kong (most 
negative with dark red color but hardly visible due to its 
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small size), Chinese Taipei, S.Korea, Vietnam. Thus the 
Google matrix approach bring a new perspective for anal¬ 
ysis of complex of economical relations between countries 
and sectors. 




panel shows the case when Be is determined by GheiRank and 
PageRank vectors; the special values are dB^s/dan = —0.0367 
for ROW group (gray) and dBwjdcjw — 0.3248 for Germany 
(magenta). Bottom panel shows the case when Be is com¬ 
puted from the Export-Import value; the special values are 
dBss/dcTii = —0.0280 fro ROW group (gray) and dBn/dan = 
0.4911 for Germany (magenta). Names of the countries can be 
found in Table and in the world map of countries [22] . 


Another results for the effects of labor cost in Germany 
and in USA are shown in Fig, and Fig.[^ In the case of 
Germany the most strong negative sensitivity is for USA, 
Russia, China for CheiRank-PageRank balance while for 
Import-Export it is Switzerland and Austria. However, 
USA and Russia are relatively weakly affected. This again 
stresses the qualitative difference between these two ap¬ 
proaches. 

The increase of USA labor cost in Fig. produces 
positive derivatives of CheiRank-PageRank balance for 
Canada and Mexico that looks reasonable from a view 
point of economy since these countries will profit from 
higher production costs in USA. In opposite, Export-Import 
gives most strong negative derivatives for Canada and 
Mexico. 

The whole matrix of labor cost derivatives dBcjdac' of 
the CheiRank-PageRank balance Be is shown in Fig. 
(numerical values of derivatives are given at [26]). Of course, 
the diagonal terms have the strongest positive derivatives. 



Fig. 20. Same as in Fig. 
labor cost c 


18 


with the derivative dBejdae' over 
= 34 of USA for year 2008. Top panel shows 


the case when Be is determined by GheiRank and PageR¬ 
ank vectors; the special values are dB^s/das^ = —0.0257 for 
ROW group (gray) and dBs^/das^ = 0.3148 for USA (ma¬ 
genta) . Bottom panel shows the case when Be is computed from 
the Export-Import value; the special values are dB^^jda^^ — 
—0.0632 for ROW group (gray) and dB^^/das^ = 0.4852 for 
USA (magenta). Names of the countries can be found in Table 
l^and in the world map of countries [22]. 


but off-diagonal terms change signs and characterize the 
sensitivity of one country to labor cost in other country. 
The vertical lines with high derivative values correspond 
to Germany (U = 11), Japan (U = 18), S.Korea (U = 19), 
USA (U = 34), China (U = 37), Russia (U = 41). The 
rest of the world (ROW) group also have a visible effect 
of other countries (U = 58). Thus is it desirable to obtain 
individual OECD data for countries of the ROW group. 


In Fig. [^ we considered the effects of the labor cost in 
various countries. We can also see the effect of price vari¬ 
ation 6s' in a given sector s' on the CheiRank-PageRank 
balance Be of country c. This sensitivity is given by the 
rectangular matrix of derivatives dBejdds' shown in Fig. 22 
(numerical data are given at [26]). The strongest pos¬ 
itive derivatives (blue squares) are for s' = 2,c = 50 
(mining and Saudi Arabia), s' = 23,c = 44 (motors and 
Hong Kong), s' = 27, c = 20 (finance and Luxembourg). 
The strongest negative derivatives (red squares) are for 
s' = 2,c = 3 (mining and Belgium), s' = 2,c = 42 (min¬ 
ing and Singapore which economy is very sensitive to min¬ 
ing products), s' = 7, c = 11 (petroleum and Germany), 
s' = 7, c = 18 (petroleum and Japan), s' = 7, c = 37 
(petroleum and China), s' = 11, c = 34 (manufacture of 
basic metals and USA), s' = ll,c = 42 (manufacture of 
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10 20 30 40 50 

C’ 


Fig. 21. Global view of the effect of labor cost variation in 
country d on country c in 2008. Matrix elements dBddcfc' 
are given in colors shown by the truncated color scale; matrix 
elements above the scale (diagonal terms) are shown in the top 
inset where y — dBcldac>- In the matrix of derivatives shown 
by color, a:-axis shows the index c of country where a labor 
cost variation (Jc> takes place and y-axis shows the country c 
affected by the change. Here Be is computed from GheiRank 
and PageRank probabilities. Gountry identification numbers 
c = 1,..., 58 are given in Table 


basic metals and Singapore). All these results are in agree¬ 
ment with the economic realities of sensitivity of the above 
countries to given activity sectors. This shows the strength 
of the Google matrix approach to analysis of WNEA. 


3.8 World transformation matrix of activity sectors 


From the obtained Google matrices G, G* of WNEA we 
can analyze the transformation of the activity sectors by 
the world economy. For this analysis we compute the trans¬ 
fer matrix 

T = (1 - r,)(l - , (15) 

where 7 ^ is a numerical constant. Our study show that 
as in the case of damping factor a the results are robust 
to variations of r] in the range 0.5 < r] < 0.9 and thus 
in the following we present the results for r] = 0.7. We 
note that a similar construction for ImpactRank has been 
used for Wikipedia networks m and the C.elegans neural 
network [ST]. In a certain sense (15) can be considered 


as a scattering matrix of particles entering in a system 
by G term and then going out by the expansion term 
1 + 77 G* + (? 7 G*)^.... = 1/(1 — ^G*). In this approach 77 
describes a relaxation rate in the system. We note that T 
belongs to the Google matrix class. 



1 0.08 


1 0.06 


0O.O4 


0.02 


I 


-0.02 


Fig. 22. Global view of the effect of sector s' price variation on 
balance of country c in 2008. Golors are proportional to matrix 
elements dBddds ', x-axis shows the sector index s' (sectors are 
given in Table and y-axis gives the country index c affected 
by the change (countries are given in Table [^. Here Be is 
computed from GheiRank and PageRank probabilities. 


From the global matrix T of size N we obtain the 
reduced matrix Rss'{c) of size Ns describing the trans¬ 
formation for activity sectors for a country c. We have 
Rss'{c') = ,c,c' where c' is a target country we are 

interested in. The matrices Rss'{c') giving the transfor¬ 
mation of sector s' to all other sectors s for c' of China, 
USA, Germany are given in [26]. The reduced transforma¬ 
tion matrix for the whole world is obtained by averaging 
over countries with Rgs' = /^d (see Fig. 23). 

The results of Fig. [^ show a few characteristic features: 
the reduced transfer matrix has a strong diagonal element 
(this is because each product is strong projection on it¬ 
self), there are characteristic horizontal lines correspond¬ 
ing to important sectors (e.g. s = 2, 7,11, 25). 


By considering a transformation of a given sector to all 
other sectors for a given country. For s' = 2 (mining) we 
present the resulting transformed vector v{s) in Fig. 
for France, Germany, Switzerland and USA. The global 
profiles are similar but there are significant enhancement 
for Germany at sector s = 1 (petroleum) and for Switzer¬ 
land at sector 5 = 20 (manufacturing and recycling). For 
comparison we show the results of transformation of in¬ 
put/output matrix M of 0 - The comparison shows a 
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10 20 30 


Fig. 23. Image of the average reduced transfer matrix Rg^s' of 
sectors to sectors for for the whole world (averaged over coun¬ 
tries) for year 2008. Here x-axis represents the initial sector s' 
and y-axis represents the hnal sectors s into which s' is trans¬ 
formed. The sector numbering is given in Table Colors are 
proportional to matrix elements and r] = 0.7. 



s 


Fig. 24. Top panel: Examples of prohle v{S) for transforma¬ 
tion vector from the reduced transfer matrix for several coun¬ 
tries in 2008. Here the initial sector is s = 2 (mining) while 
the transformed vector x(s) is formed by the matrix defined 
in Eig. the countries are Erance (blue), Germany (red), 
Switzerland (green) and USA (black). Bottom panel: Eor com¬ 
parison, we show here the same as top panel but instead of T, R 
matrices we use the input/output matrix M with normalized 
columns (dangling nodes are not replaced here, transitions in¬ 
side one country are taken to be zero); a column s' of such a 
matrix for country c' is given by ,cc'] here the same 

countries are shown by same colors as in top panel.. 


drastic difference between two approaches which we at¬ 
tribute to the fact that M does not take into account the 
multiple network transitions. 

The transformation for the sector s' = 34 are shown in 
Fig. 25 for Cyprus (blue), Singapore (red), Luxembourg 
(green) and Malta (black). We see that for Luxembourg 
there is a strong transformation of s' = 20 to 5 = 6 
(publishing). At the same time the global profile, being 
different from the case of Fig. 24 with s' = 2, has simi¬ 
lar features for different countries. The comparison with 
the transformation results from value exchange matrix 
^ss',cc' are again very different as in the case of Fig. [24 

The obtained results for the activity sector transto- 
mation by the WNEA open new possibilities for analysis 
of interactions between the world economic activities. The 
Google matrix approach provides new type of results being 
very different from usual Input/Output matrix approach. 
This is related to the fact that the transformation matrix 
(14) takes into account summation over various cycles over 
the network. 



(red), Luxembourg (green) and Malta (black). 


4 Discussion 

In this work we have developed the Google matrix anal¬ 
ysis of the world network of economic activities from the 
OECD-WTO TiVA database. The PageRank and CheiRank 
probabilities allowed to obtain ranking of world countries 
independently of their richness being mainly determined 
by the efficiency of their economic relations. The devel¬ 
oped approach demonstrated the asymmetry in the eco¬ 
nomic activity sectors some of which are export oriented 


and others are import oriented. We also showed that the 
eigenstates of the WNEA Google matrix select specific 
quasi-isolated communities oriented to specific activity sec¬ 
tors. The CheiRank-PageRank balance Be allows to deter¬ 
mine economically rising countries with robust network 
of economic relations. The sensitivity of this Be to price 
variations and labor cost in various countries determines 
the hidden relations between world economies being not 
visible via usual Export-Import exchange analysis. The 
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Google matrix analysis determines also the transforma¬ 
tion features of world activity sectors. 

The comparison with the multiproduct world trade 
network from UN COMTRADE shows certain similari¬ 
ties between the two networks of WNEA and WTN. At 
the same time the WNEA data provides new elements 
for interactions of activity sectors while there are no di¬ 
rect interactions of products in COMTRADE database. 
From this viewpoint the OECD-WTO data captures the 
economic reality on a deeper level. But at the same time 
the OECD-WTO network is less developed compared to 
COMTRADE (less countries, years, sectors). Thus it is 
highly desirable to extend the OECD-WTO database. 

We think that the Google matrix analysis developed 
here and in |13ll4j captures better the new reality of mul¬ 
tifunctional directed tensor interactions and that the uni¬ 
versal features of this approach can be also extended to 
multifunctional financial network flows which now attract 
an active interest of researchers USES]. Unfortunately, the 
data on financial flows have much less accessibility com¬ 
pared to the networks discussed here. 

We point that recently some of the matrix methods, 
developed in physics community, started to find active ap¬ 
plication for economy systems (see e.g. mm)- However, 
usually for physicists these matrices have been from the 
unitary or Hermitian ensembles, where the Random Ma¬ 
trix Theory allowed to obtained certain universal results. 
Here, we show that the directed networks and tensors ap¬ 
pearing in the interacting economy systems are described 
by the matrices of Perron-Frobenius operators which had 
not been studied much in physics. Thus the new field of 
research is now opened for physicists, mathematicians and 
computer scientists with application to complex interact¬ 
ing economy systems. 
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country name 

country code 

country flag 


country name 

country code 

country flag 

1 

Australia 

AUS 


30 

Sweden 

SWE 


2 

Austria 

AUT 


31 

Switzerland 

CHE 

o 

3 

Belgium 

BEL 

U 

32 

Turkey 

TUR 


4 

Canada 

CAN 


33 

United Kingdom 

GBR 


5 

Chile 

CHL 


34 

United States 

USA 

# 

6 

Czech Republic 

CZE 


35 

Argentina 

ARG 

/< 

o ^ 

7 

Denmark 

DNK 


36 

Brazil 

BRA 


8 

Estonia 

EST 


37 

China 

GHN 


9 

Finland 

FIN 


38 

Chinese Taipei 

TWN 

• 

10 

France 

FRA 

o 

39 

India 

IND 

ASk. 

11 

Germany 

DEU 

• 

40 

Indonesia 

IDN 


12 

Greece 

GRC 


41 

Russia 

RUS 

IP 









13 

Hungary 

HUN 


42 

Singapore 

SGP 


14 

Iceland 

ISL 


43 

South Africa 

ZAP 


15 

Ireland 

IRL 

o 

44 

Hong Kong 

HKG 

Q 

16 

Israel 

ISR 

✓ X 

45 

Malaysia 

MYS 


17 

Italy 

ITA 

o 

46 

Phillippines 

PHL 

d 

18 

Japan 

JPN 


47 

Thailand 

THA 


19 

Korea 

KOR 


48 

Romania 

ROU 

U 

20 

Luxembourg 

LUX 

s 

49 

Vietnam 

VNM 

o 

21 

Mexico 

MEX 

u 

50 

Saudi Arabia 

SAU 

Q 









22 

Netherlands 

NLD 


51 

Brunei Darussalam 

BRN 


23 

New Zealand 

NZL 

• 

52 

Bulgaria 

BGR 

w 

24 

Norway 

NOR 


53 

Cyprus 

CYP 










25 

Poland 

POL 


54 

Latvia 

LVA 


26 

Portugal 

PRT 


55 

Lithuania 

LTU 

• 

27 

Slovak Republic 

SVK 


56 

Malta 

MLT 


28 

Slovenia 

SVN 


57 

Cambodia 

KHM 

Q 

29 

Spain 

ESP 

e 

58 

Rest of the World 

ROW 

9 


Table 1. List of Nc = 58 countries (with rest of the world ROW) with country name, code and flag. 
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OECD ICIO Category 

ISIC Rev. 3 correspondence 

1 

C01T05 AGR 

01 - Agriculture, hunting and related service activities 

02 - Forestry, logging and related service activities 

05 - Fishing, operation of fish hatcheries and fish farms; service activities incidental to fishing 

2 

C10T14 MIN 

10 - Mining of coal and lignite; extraction of peat 

11 - Extraction of crude petroleum and natural gas; service activities incidental to oil and gas extraction excluding surveying 

12 - Mining of uranium and thorium ores 

13 - Mining of metal ores 

14 - Other mining and quarrying 

3 

C15T16 FOD 

15 - Manufacture of food products and beverages 

16 - Manufacture of tobacco products 

4 

C17T19 TEX 

17 - Manufacture of textiles 

18 - Manufacture of wearing apparel; dressing and dyeing of fur 

19 - Tanning and dressing of leather; manufacture of luggage, handbags, saddlery, harness and footwear 

5 

C20 WOD 

20 - Manufacture of wood and of products of wood and cork, except furniture; 

Manufacture of articles of straw and plaiting materials 

6 

C21T22 PAP 

21 - Manufacture of paper and paper products 

22 - Publishing, printing and reproduction of recorded media 

7 

C23 PET 

23 - Manufacture of coke, refined petroleum products and nuclear fuel 

8 

C24 CHM 

24 - Manufacture of chemicals and chemical products 

9 

C25 RBP 

25 - Manufacture of rubber and plastics products 

10 

C26 NMM 

26 - Manufacture of other non-metallic mineral products 

11 

C27 MET 

27 - Manufacture of basic metals 

12 

C28 FBM 

28 - Manufacture of fabricated metal products, except machinery and equipment 

13 

C29 MEQ 

29 - Manufacture of machinery and equipment n.e.c. 

14 

C30 ITQ 

30 - Manufacture of office, accounting and computing machinery 

15 

C31 ELQ 

31 - Manufacture of electrical machinery and apparatus n.e.c. 

16 

C32 CMQ 

32 - Manufacture of radio, television and communication equipment and apparatus 

17 

C33 SCQ 

33 - Manufacture of medical, precision and optical instruments, watches and clocks 

18 

C34 MTR 

34 - Manufacture of motor vehicles, trailers and semi-trailers 

19 

C35 TRQ 

35 - Manufacture of other transport equipment 

20 

C36T37 OTM 

36 - Manufacture of furniture; manufacturing n.e.c. 

37 - Recycling 

21 

C40T41 EGW 

40 - Electricity, gas, steam and hot water supply 

41 - Collection, purification and distribution of water 

22 

C45 CON 

45 - Construction 

23 

C50T52 WRT 

50 - Sale, maintenance and repair of motor vehicles and motorcycles; retail sale of automotive fuel 

51 - Wholesale trade and commission trade, except of motor vehicles and motorcycles 

52 - Retail trade, except of motor vehicles and motorcycles; repair of personal and household goods 

24 

C55 HTR 

55 - Hotels and restaurants 

25 

C60T63 TRN 

60 - Land transport; transport via pipelines 

61 - Water transport 

62 - Air transport 

63 - Supporting and auxiliary transport activities; activities of travel agencies 

26 

C64 PTL 

64 - Post and telecommunications 

27 

C65T67 FIN 

65 - Financial intermediation, except insurance and pension funding 

66 - Insurance and pension funding, except compulsory social security 

67 - Activities auxiliary to financial intermediation 

28 

C70 REA 

70 - Real estate activities 

29 

C71 RMQ 

71 - Renting of machinery and equipment without operator and of personal and household goods 

30 

C72 ITS 

72 - Computer and related activities 

31 

C73 RDS 

73 - Research and development 

32 

C74 BZS 

74 - Other business activities 

33 

C75 GOV 

75 - Public administration and defense; compulsory social security 

34 

C80 EDU 

80 - Education 

35 

C85 HTH 

85 - Health and social work 

36 

C90T93 OTS 

90 - Sewage and refuse disposal, sanitation and similar activities 

91 - Activities of membership organizations n.e.c. 

92 - Recreational, cultural and sporting activities 

93 - Other service activities 

37 

C95 PVH 

95 - Private households with employed persons 


Table 2. List of sectors considered by Input/Output matrices from OECD database, their correspondence to the ISIC classifi¬ 
cation is also given. 




20 


V.Kandiah, H.Escaith and D.L.Shepelyansky: Google matrix of the world network of economic activities 


Sector 

K (1995) 

% vol (1995) 

k* (1995) 

% vol (1995) 

k (2008) 

% vol (2008) 

k* (2008) 

% vol (2008) 

1 

19 

2.2979 

16 

2.9763 

20 

1.9532 

16 

2.0902 

2 

27 

1.2993 

2 

8.6183 

24 

1.5245 

1 

15.8784 

3 

3 

6.0117 

12 

3.3271 

11 

3.9327 

17 

1.9835 

4 

10 

3.9579 

14 

3.0831 

17 

2.0934 

19 

1.8634 

5 

30 

1.108 

20 

1.9037 

33 

0.60075 

22 

1.3001 

6 

11 

3.5687 

6 

4.2128 

18 

2.0608 

14 

2.3736 

7 

4 

5.9126 

19 

2.2783 

1 

11.589 

4 

6.34 

8 

2 

6.251 

1 

10.6954 

3 

6.0558 

2 

9.1103 

9 

17 

2.4035 

15 

3.0546 

19 

1.9785 

13 

2.5549 

10 

28 

1.2 

21 

1.8337 

29 

1.0389 

21 

1.3177 

11 

8 

4.4393 

3 

8.0658 

4 

5.4907 

3 

8.3184 

12 

20 

2.2646 

17 

2.7194 

23 

1.6212 

15 

2.2182 

13 

9 

4.0642 

8 

4.0365 

9 

4.0117 

9 

4.0597 

14 

12 

3.3353 

13 

3.158 

8 

4.0642 

6 

5.0066 

15 

18 

2.3789 

9 

4.0148 

25 

1.456 

18 

1.8673 

16 

15 

2.7053 

10 

3.8054 

14 

2.7844 

11 

3.6339 

17 

31 

1.0034 

23 

1.1434 

34 

0.31041 

29 

0.40161 

18 

7 

5.2722 

7 

4.1643 

6 

5.1478 

10 

3.9907 

19 

26 

1.3665 

22 

1.7813 

26 

1.3028 

23 

1.2752 

20 

24 

1.6331 

27 

0.67546 

22 

1.6652 

20 

1.3858 

21 

21 

2.1673 

30 

0.34377 

10 

3.946 

30 

0.39969 

22 

1 

6.538 

32 

0.22022 

2 

6.8692 

32 

0.15209 

23 

5 

5.8472 

4 

7.9296 

7 

4.6893 

8 

4.6745 

24 

25 

1.5283 

29 

0.37682 

27 

1.2377 

27 

0.62202 

25 

6 

5.8385 

5 

6.5023 

5 

5.2454 

5 

5.8065 

26 

29 

1.1862 

26 

0.6839 

28 

1.2179 

26 

0.62929 

27 

13 

2.7584 

18 

2.3006 

15 

2.5623 

12 

3.3487 

28 

33 

0.70446 

24 

0.93849 

31 

0.84772 

33 

0.105 

29 

36 

0.16329 

33 

0.18955 

36 

0.21276 

24 

0.81082 

30 

34 

0.53799 

28 

0.39581 

32 

0.67481 

28 

0.61668 

31 

35 

0.36919 

31 

0.33351 

35 

0.24684 

31 

0.24177 

32 

16 

2.618 

11 

3.372 

13 

3.0455 

7 

4.7163 

33 

14 

2.7071 

34 

0.064931 

12 

3.3939 

35 

0.06377 

34 

32 

0.89993 

36 

0.0416 

30 

1.036 

34 

0.09439 

35 

22 

1.8912 

35 

0.045551 

16 

2.2601 

36 

0.025979 

36 

23 

1.7326 

25 

0.7136 

21 

1.8131 

25 

0.72283 

37 

37 

0.03899 

37 

0 

37 

0.019524 

37 

0 


Table 3. First column gives the sectors from OECD database, for each of them the following columns give the Import Rank K 
with the sector fraction in global trade value and ExportRank X* with sector fraction in global trade value. Data are shown 
for 1995 and 2008. 
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K 

K- 

K2 

k 


1 

DEU C34 MTR 

ROW C10T14 MIN 

DEU C24 CHM 

USA C23 PET 

ROW C10T14 MIN 

2 

USA C75 GOV 

RUS C10T14 MIN 

USA C65T67 EIN 

JPN C23 PET 

SAU C10T14 MIN 

3 

ROW C75 GOV 

SAU C10T14 MIN 

DEU C29 MEQ 

USA C75 GOV 

RUS C10T14 MIN 

4 

SAU C85 HTH 

USA C24 CHM 

DEU C34 MTR 

ROW C45 CON 

USA C24 CHM 

5 

GBR C85 HTH 

DEU C24 CHM 

DEU C27 MET 

CHN C32 CMQ 

CAN C10T14 MIN 

6 

USA C34 MTR 

DEU C27 MET 

USA C74 BZS 

CHN C27 MET 

DEU C24 CHM 

7 

ROW C45 CON 

NOR C10T14 MIN 

DEU C50T52 WRT 

USA C45 CON 

NOR C10T14 MIN 

8 

ROW C15T16 EOD 

RUS C27 MET 

USA C24 CHM 

DEU C34 MTR 

AUS C10T14 MIN 

9 

USA C15T16 EOD 

USA C50T52 WRT 

DNK C60T63 TRN 

KOR C23 PET 

CHN C30 ITQ 

10 

RUS C50T52 WRT 

DEU C29 MEQ 

GBR C74 BZS 

DEU C23 PET 

USA C30 ITQ 

11 

USA C45 CON 

USA C74 BZS 

JPN C34 MTR 

JPN C40T41 EGW 

JPN C30 ITQ 

12 

USA C85 HTH 

CHN C27 MET 

GBR C65T67 EIN 

ROW C75 GOV 

DEU C29 MEQ 

13 

DEU C15T16 EOD 

USA C60T63 TRN 

CHN C32 CMQ 

CHN C24 CHM 

DEU C34 MTR 

14 

ROW C60T63 TRN 

GBR C65T67 EIN 

CHN C24 CHM 

USA C34 MTR 

KOR C32 CMQ 

15 

USA C65T67 EIN 

USA C23 PET 

DEU C60T63 TRN 

USA C24 CHM 

USA C23 PET 

16 

GBR C50T52 WRT 

GBR C74 BZS 

ERA C50T52 WRT 

CHN C30 ITQ 

USA C74 BZS 

17 

DEU C24 CHM 

USA C65T67 EIN 

USA C50T52 WRT 

CHN C23 PET 

TWN C32 CMQ 

18 

DEU C29 MEQ 

CHN C30 ITQ 

CHN C50T52 WRT 

ROW C60T63 TRN 

CHN C27 MET 

19 

DEU C50T52 WRT 

DEU C34 MTR 

CHN C29 MEQ 

CHN C29 MEQ 

DEU C27 MET 

20 

DEU C27 MET 

USA C30 ITQ 

ROW C60T63 TRN 

DEU C29 MEQ 

GBR C74 BZS 


Table 4. Top 20 ranks for global PageRank CheiRankK*, 2DRank i^ 2 , ImportRank K and ExportRank K* for the year 
2008. 


Ki 

IV’il 

node 

iV’il 

node 

iV’ii 

node 

IV'il 

node 

1 

0.037606 

ROW C17T19 TEX 

0.050431 

ARC C34 MTR 

0.054681 

CHN C32 CMQ 

0.052248 

RUS C10T14 MIN 

2 

0.025695 

CHN C17T19 TEX 

0.049991 

BRA C34 MTR 

0.053306 

KOR C32 CMQ 

0.03948 

SAU C10T14 MIN 

3 

0.021618 

ITA C17T19 TEX 

0.029753 

JPN C34 MTR 

0.053253 

TWN C32 CMQ 

0.026187 

ROW C10T14 MIN 

4 

0.017075 

USA C17T19 TEX 

0.026592 

DEU C34 MTR 

0.027361 

SGP C32 CMQ 

0.022125 

NOR C10T14 MIN 

5 

0.016216 

CHN C32 CMQ 

0.018372 

THA C34 MTR 

0.025189 

MYS C32 CMQ 

0.019764 

USA C71 RMQ 

6 

0.013003 

CHN C30 ITQ 

0.01531 

IDN C34 MTR 

0.018824 

USA C30 ITQ 

0.013899 

USA C50T52 WRT 

7 

0.010963 

ERA C17T19 TEX 

0.0093875 

ROW C21T22 PAP 

0.016965 

PHL C32 CMQ 

0.011638 

ROW C29 MEQ 

8 

0.010175 

TUR C17T19 TEX 

0.0093382 

DEU C15T16 EOD 

0.01534 

JPN C30 ITQ 

0.010871 

RUS C27 MET 

9 

0.010161 

USA C75 GOV 

0.0090288 

USA C15T16 EOD 

0.014664 

GBR C65T67 EIN 

0.0082943 

DEU C29 MEQ 

10 

0.0099839 

USA C65T67 EIN 

0.0086552 

USA C75 GOV 

0.013713 

CHN C30 ITQ 

0.0082905 

RUS C23 PET 


Table 5. Top 10 values of 4 different eigenvectors from Figlsj : „ _ 

A = 0.4993 (red), A = 0.3746 + 0.0126z (green), A = 0.6256 (blue) and A = —0 
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The 


corresponding eigenvalues from left to right are 
.0001 + 0.1687z (magenta). 




























